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WHAT IS CLAIMED IS: 

1. A method for numerical analysis of a flow field 
of incompressible viscous fluid/ directly using V-CAD data, 
5 comprising: 

a dividing step (A) of dividing external data into a 
plurality of cells (13) having boundaries orthogonal to 
each other, the external data including boundary data of an 
object which contacts incompressible viscous fluid; 
10 a cell classifying step (B) of classifying the 

divided cells into an internal cell (13a) positioned inside 
or outside the object and a boundary ceil (13b) including 
the boundary data; 

a cut point determining step (C) of determining cut 
15 points in ridges of the boundary cell on the basis of the 
boundary data; 

a boundary face determining step (D) of determining a 
polygon connecting the cut points to be cell internal data 
for the boundary face; and 
20 . a analyzing step (E) of applying a cut cell finite 

volume method combined with a VOF method to a boundary of a 
flow field to analyze the flow field. 



2. A method for numerical analysis of a flow field 
25 of incompressible viscous fluid, according to claim 1, the 
analyzing step (E) comprising: 
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applying a two-dimensional QUICK interpolation scheme 
to a convection term for space integral; 

applying central difference having precision of a 
degree of a second order to a diffusion term; 
5 combining the convection term and the diffusion term, 

and applying Adams-Bashf orth method having precision of a 
degree of a second order to the combined convection term 
and diffusion term for time marching; and 

applying an Euler implicit method having precision of 
10 a degree of a first order to a pressure gradient term for 
time marching. 

3. A method for numerical analysis of a flow field 
"o"f incompressible viscous fluid, according to claim 2, 

15 wherein for a two-dimensional boundary cell, a governing 
equation in the finite volume method is expressed by a 
governing equation (7) of Formula 1, 
[Formula 1] 

jj^V = - Ijdiviu ® iiW - jjdiv{pI)dV + j\div(grad{u))ciV ( 7 ) 

20 

4. A method for numerical analysis of a flow field 
of incompressible viscous fluid, according to claim 3, 
wherein the convection term, the pressure gradient term and 
the diffusion term in the governing equation of the finite 

25 volume method are expressed by the equations (8),. (9) and 
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(10) of Formula 2, respectively, 
[Formula 2] 

convection term: 

=[w.;":x-5,-.,"i:>,-,,,) 

-K Ay(5, .v<'> M, ,. - B. , ,.v'" , ) 

-^dLxiAf jv\^]v. J - 4 j-i^/,y-i^/,;-i )]7 \only no - o« wall 
pressure gradient term: 

Jpiv(/>/)t/K=cg(p/) H</S== X /'^^'^^^n. 
diffusion term: 

5. A method for numerical analysis of a flow field 
of incompressible viscous fluid, according to claim 3, 
wherein when a no-slip boundary condition is used for a 
solid boundary, integral is performed on the solid boundary 
with the convection term being zero, a value of a middle 
point P of a cut line segment being used as an average for 
the pressure gradient term and the diffusion term, and a 
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space integral is performed with areas fractions being 
applied to all of the terms. 

6. A method for numerical analysis of a flow field 
5 of incompressible viscous fluid, according to claim 3, 

wherein the boundary cell having the parameter smaller than 
a threshold value of VOF=0.01 is regarded as a complete 
solid, 

for the boundary cell having the parameter larger 
10 than the threshold value, a definition point for the 

parameter calculated in a cut cell is set at a center of 
the boundary cell, 

and a definition point for a parameter in a ridge is 
set at a center of a cell ridge, and a parameter at a 
15 middle point of a line segment 4 is calculated by a linear 
interpolation . 

7. A method for numerical analysis of a flow field 
of incompressible viscous fluid, according to claim 3, 

20 wherein a drag force (in a flow direction) and a lift force 
(in a direction vertical to the flow) acting on the object 
are expressed by equations (12) and (13) of Formula 3, 



25 
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[Formula 3] 



drag force: 



y 



lift force: 



Cartesian 



y 

y2 



V^y^ l/.(.v) U(.v) J(CT^ l^^^.j -O-^ Wu>>&|on/, 



Cartesian 



(12) 



(13) 



8. A method for numerical analysis of a flow field 
of incompressible viscous fluid, according to claim. 3, 
wherein in fluid-structure interaction analysis 
accompanying a moving boundary, a fluid system and a 
structure system are separately analyzed each predetermined 
time interval, and boundary conditions for the fluid system 
and the structure system are explicitly used. 



9* A method for numerical analysis of a flow field 
of incompressible viscous fluid, according to claim 8, 
wherein in analysis on a forcibly vibrated circular 
cylinder, the circular cylinder is set as one-mass-point 
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and one-degree-of-freedom system such that the circular 
cylinder is a solid structure elastically supported and 
vibrating in a direction vertical to the flow, 

and Y-direction displacement of a center of the 
circular cylinder is given by the equation (17), and a 
velocity boundary condition in the Y direction for a 
surface of the circular cylinder is given by the equation 
(18), 

y=Asin (2 7C fct) ...(17) 

Vw=A2 7i:fcCOs(2 7cfct) ...(18). 

10. A method for numerical analysis of a flow field 
Of incompressible viscous fluid, according to claim 9, 
wherein movement velocity of the vibrating circular 
cylinder obtained by the equation (18) is changed to be 
given each calculation time step for the velocity boundary 
condition on the flow field. 

11. A method for numerical analysis of a flow field 
of incompressible viscous fluid, according to claim 8, 
wherein in analysis on self-induced vibration due to an 
vortex shedding from the circular cylinder, a vibration 
equation having a dimension is expressed by the equation 
(19) or (20), using one-mass-point and one-degree-of- 
freedom dumper /spring model, 
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m 



dt^ dt 1 



(19) 




(20) 



12. A method for numerical analysis of a flow field 
of incompressible viscous fluid, according to claim 11, 
wherein the movement velocity of the vibrating circular 
cylinder calculated by the equation (20) is changed to be 
given each calculation time step for the velocity boundary 
condition on the flow field. 

13. A method for numerical analysis of a flow field 
of incompressible viscous fluid, according to claim 11, 
wherein initial displacement and initial velocity of the 
circular cylinder are set to be zero, the lift force is 
explicitly given by using a current value, and the 
vibration equation is integral by the Newmark' s /3 method 
to obtain vibration displacement and vibration velocity of 
the circular cylinder. 

14. A device for numerical analysis of a flow field 
of incompressible viscous fluid, directly using V-CAD data, 
comprising: 

an input device (2) for inputting external data 
including boundary data of an object (1) that contacts 
incompressible viscous fluid; 
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an external storage device (3) for storing 
substantial data of shape data and physical property data 
integrated into each other, and a storage operational 
program for the substantial data; 
5 an internal storage device (4) and central processing 

device (5) for executing the storage operational program; 
and 

an output device (6) for outputting a result of the 
execution; 

10 wherein the device divides the external data into a 

plurality of cells (13) having boundaries orthogonal to 
each other, classifies the divided cells into an internal 
cell (13a) positioned inside or outside the object and a 
boundary cell (13b) including the boundary data, determines 

15 cut points in ridges of the boundary cell on the basis of 
the boundary data, determines a polygon connecting the cut 
points to be cell internal data for the boundary face, and 
applies a cut cell finite volume method combined with a VOF 
method to a boundary of a flow field to analyze the flow 

20 field. 

15. A program for numerical analysis of a flow field 
of incompressible viscous fluid, directly using V-CAD data, 
causing a computer to perform: 
25 a dividing step (A) of dividing external data into a ' 

plurality of cells (13) having boundaries orthogonal to 
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each other, the external data including boundary data of an 
object which contacts incompressible viscous fluid; 

a cell classifying step (B) of classifying the 
divided cells into an internal cell (13a) positioned inside 
or outside the object and a boundary cell (13b) including 
the boundary data; 

a cut point determining step (C) of determining cut 
points in ridges of the boundary cell on the basis of the 
boundary data; 

a boundary face determining step (D) of determining a 
polygon connecting the cut points to be cell internal data 
for the boundary face; and 

a analyzing step (E) of applying a cut cell finite 
volume method combined with a VOF method to a boundary of a 
flow field to analyze the flow field. 



